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A ubiquitous observation in crowded cell membranes is that molecular transport does not follow Fickian diffusion but exhibits subdiffusion. 

The microscopic origin of such a behaviour is not understood and highly debated. Here we discuss the spatio-temporal dynamics for two 
(~~) models of subdiffusion: fractional Brownian motion and hindered motion due to immobile obstacles. We show that the different microscopic 
y— H mechanisms can be distinguished using fluorescence correlation spectroscopy (FCS) by systematic variation of the confocal detection area. 

We provide a theoretical framework for space-resolved FCS by generalising FCS theory beyond the common assumption of spatially 

Gaussian transport. We derive a master formula for the FCS autocorrelation function, from which it is evident that the beam waist of an 
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FCS experiment is a similarly important parameter as the wavenumber of scattering experiments. These results lead to scaling properties 
of the FCS correlation for both models, which are tested by in silico experiments. Further, our scaling prediction is compatible with the 
FCS half-value times reported by Wawrezinieck et al. [Biophys. J. 89, 4029 (2005)] for in vivo experiments on a transmembrane protein. 

Key words: single-molecule techniques, diffusion of macromolecules, random processes, molecular crowding 
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Introduction 

Measurement of molecular transport at the subcellular level can 
provide important information on both physiological mechanisms 
and physical interactions that drive and constrain biochemical pro- 
cesses. The obstructed motion of biomolecules in living cells dis- 
plays anomalous transport including subdiffusion, which was es- 
tablished in the past decade by numerous experiments applying 
techniques with labeled particles. Nevertheless, the interpretation 
of the collected data remains often controversial and the origin 
of the subdiffusive behaviour is highly debated. 1-6 Crowded en- 
vironments like cellular membranes contain structures on many 
length scales, and further progress depends on experimental tech- 
niques that resolve transport on these different scales. Such spatio- 
temporal information is needed to test and refine models of anoma- 
lous transport. 

One widespread technique for the investigation of molecular 
transport is fluorescence correlation spectroscopy (FCS), which 
follows the motion of fluorescently labeled molecules with high 
temporal resolution. 7,8 This mesoscopic, local method consists of 
collecting the fluorescent light from a steadily illuminated volume 
or area and autocorrelating its intensity fluctuations. An important 
parameter of FCS measurements is the beam waist of the illumina- 
tion laser. While experimental setups in the past were constrained 
to a fixed value, recent technological advancements allow large 
variations of the confocal detection area to gather spatial informa- 
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tion.- — By z-scan FCS^i^ or beam expanders, 11 focal radii can 
be varied between 200 nm and 500 nm, but measurements at the 
nanoscale became possible by introducing nano-apertures (75 nm 
to 250nm). 12 Only lately, a far-field optical nanoscopic method 
named stimulated emission depletion (STED) fluorescence corre- 
lation microscopy was developed to beat the diffraction limit,— al- 
lowing focal radii to span almost a decade down to 15 nm (Ref. 
0). 

The subdiffusive motion of macromolecules in crowded cells 
and membranes was studied extensively by FCS experiments^— 
and was complemented in real-space by single-particle track- 
ing,— — For Fickian diffusion, the mean-square displacement 
grows linearly in time, Sr 2 (i) = 4Dt in two dimensions with dif- 
fusion constant D. Then, the decay of the FCS autocorrelation 
function obeys 

1 1 

G(t) = - - , — , (1) 



N 1 + t/T D 

where r D = w 2 /4D denotes the dwell time and N the average num- 
ber of labeled molecules in the illuminated areaQ These equations 
are no longer valid for anomalous transport. Introducing the walk 
dimension rf w , subdiffusion is characterised by 6r 2 (t) ~ r 2/rfw , and 
FCS experiments are often rationalised by 



G(t)=l- 



1 



N 1 + (Tf)° (2) 
upon fitting N, T, and a. It is usually and tacitly anticipated that 
both exponents coincide, a = 2/d w . 

Here, we provide a theoretical framework for space-resolved 
FCS. Relating the FCS function G(f) to the intermediate scattering 



' We restrict the discussion to two dimensional systems relevant for membranes, 
where focus distortions are negligible; in three dimensions, the asphericity of the 
illumination volume renders the formulae more cumbersome. We also ignore ef- 
fects due to the photophysics of the dye molecules, which are relevant at very short 
time scales only. This does not effect the generality of our discussion nor any of 
our conclusions. 
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function, we generalise the conventionally used fit models and con- 
nect FCS to time-resolved scattering techniques. If the beam waist 
is considered an adjustable experimental parameter similar to the 
scattering angle, FCS is turned into a valuable tool for the investi- 
gation of complex and in particular anomalous transport. The new 
approach greatly facilitates in silico experiments: for two models 
of subdiffusion, we show how spatio-temporal information on the 
tracer dynamics can be obtained and used to distinguish different 
mechanisms as the origin of anomalous transport. 

Theory 

Generalised FCS theory. Let us briefly revisit the theory under- 
lying the FCS technique;-*— we specialise to two dimensions for 
simplicity. The detection area is illuminated by a laser beam with 
intensity profile W(r). The fluorescent light depends on the fluc- 
tuating, local concentration c(r, /) of labeled molecules in the laser 
focus. Thus, the intensity collected at the detector is a spatially 
weighted average, 1(f) cc f d 2 r W(r) c(r, t). The output of the FCS 
experiment is the time-autocorrelation function of the intensity 
fluctuation 51(f) = 1(f)- (I) around the mean intensity. It is conven- 
tionally normalised as G(f) = (61(f) 61(0)) / (I) 2 ; proper normalisa- 
tion would be achieved by multiplication with N = (I) 2 / (<W 2 )- 
Introducing spatial Fourier transforms, one arrives at the represen- 
tation 

1 jd 2 q \W( q )\ 2 S( q ,t) 

G(t; w) = , (3) 

N f^q\W(q)\ 2 S(q,t = 0) 

where 5(q, f) = jd 2 r exp(iq • r) (6c(r, f) 6c(0, 0)> is known as 
the intermediate scattering function and W(q) denotes the Fourier 
transform of the intensity profile W(r). 

An conventional laser emits a Gaussian beam profile, W(r) oc 
exp (-2r 2 /w 2 ), with beam waist w, which implies a Gaussian filter 
function | W(q)| 2 oc exp (-q 2 w 2 /4j. Usually only a small fraction of 
the molecules is labeled, and then 5 (q, reduces to the incoherent 
intermediate scattering function 

5 (q, t) » F (q, t) = (exp (iq • AR(f))) . (4) 

Considering the displacements AR(?) := R(f) - R(0) after a fixed 
time lag a random variable, the incoherent scattering function can 
be interpreted as their characteristic function. For Gaussian and 
isotropic displacements, (AR(?)) = 0, only the second cumulant 
6r 2 (t) : = (|AR(0l 2 ) is non-zero. Thus F(q,f) = exp (-q 2 6r 2 (t)/4) 
for two-dimensional motion. The corresponding FCS function is 
calculated to 

G Gauss ( f ;w) = i i+(5r ! (o/w2 . (5) 

For normal diffusion, it holds F(q, t) = exp [-Dq 2 tj, and G(f) at- 
tains the simple form of Eq. l[T). For the case of subdiffusion, 
tSr 2 ^) ~ f, and Gaussian spatial displacements as in fractional 
Brownian motion (FBM), one recovers the conventional expres- 
sion, Eq. {2). 

In many complex systems, however, the (strong) assumption of 
Gaussian displacements is not valid and may only serve as an ap- 
proximation. This assumption can be tested experimentally by re- 
solving the spatial properties of the particle trajectories. An exact 



expression for the FCS function is obtained by combining equa- 
tions {3]l and l@. Evaluating the integrals over the wavenumber 
yields 

G(t; w) = N~ i (exp (-AR(r) 2 /w 2 )) , (6) 

which is a central result of our work. Let us emphasise that it 
does not require any assumptions on the dynamics; corrections may 
arise from non-dilute labeling of the molecules and from deviations 
of the Gaussian beam profile. In three-dimensional systems, one 
should further correct for anisotropies in the confocal volume. This 
expression enables new insight in the potential of the FCS tech- 
nique with consequences for the design of future FCS experiments. 
The similarity of the representation of G(t; w) in Eq. {6]l with that 
of F(k, f) in Eq. © suggests that FCS encodes important spatial 
information analogous to scattering methods like photon correla- 
tion spectroscopy or neutron spin echo. In the case of anomalous 
transport discussed below, we will use it as starting point for the 
derivation of the scaling properties of G(t; w). Equation ® shows 
that the FCS function NG(t) can be neatly interpreted as the return 
probability for a fluorescent molecule to be again (or still) in the il- 
luminated area[| As a by-product, it provides a simple description 
for the efficient evaluation of autocorrelated FCS data in computer 
simulations, circumventing the evaluation of the rapidly fluctuating 
fluorescent light intensity. 

Van Hove correlation function. The dynamics of a single la- 
beled particle is encoded in the probability distribution of the time- 
dependent displacements, P(r, t) = {6(r - R(f)>; due to rotational 
symmetry, it actually depends merely on the magnitude r = |r|. 
This function is also known as van Hove (self-)correlation func- 
tion G(r, t) in the field of liquid dynamics; 2 i to avoid confusion 
with the FCS function, we follow the notation of Ref . |U Explicit 
expressions for P(r, t) exist for many models, but for the dynamics 
on percolation clusters only conjectures of the asymptotic scaling 
behaviour are available. Let us consider a random walker on the in- 
cipient infinite percolation cluster, i.e., precisely at the percolation 
threshold. Then, the dynamics is characterised by two universal 
exponents: the fractal dimension d { and the walk dimension d w . 
Let further / and a denote the typical microscopic time and length 
scales, respectively. The van Hove function is expected to obey the 
following scaling law for r » <r and t » t o (Ref. 

Poo(r,f) = r^P^rf 11 ''*) . (7) 

The subscript oo indicates that the average is taken only for tracers 
on the infinite cluster. During a time t, the walker explores re- 
gions of linear extension of the order of R ~ t 1/<,w . The probability 
for larger excursions decreases rapidly (presumably like a stretched 
exponential), hence we assume P m (x » 1) — » rapidly. This prop- 
erty specifies the time evolution of the mean-square displacement 
and of higher moments. For the FCS measurements, however, we 
additionally need the limiting behaviour of the scaling function for 
small arguments. 



For sufficiently large time lag, the probability to find the fluorophore at a particular 
point within the confocal volume becomes independent of the position. Then, the 
FCS function can be approximated by the probability of being at or returning to 
the centre of the confocal volume after the given time multiplied by the size of the 
confocal volume. 
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Fig. 1: Simulated FCS correlation function Goa(t) on a logarithmic time 
axis for tracers on the infinite cluster at the critical obstacle density; the 
beam waist w of the laser increases from left to right. Dotted lines cor- 
respond to FBM with identical mean-square displacement, Eq. {3); the 
broken line displays the FCS function for unobstructed, normal diffusion 
(w = 450 nm). Inset: subdiffusive behaviour of the mean-square displace- 
ment for the simulated obstructed motion (solid line), again compared to 
the case of normal diffusion (broken line). 

Return probability. Integrating the van Hove function over dis- 
tances r < w with w much larger than any microscopic length 
yields the probability n(f; w) to return to the starting point of 
the random walk within a radius w after a time t. Provided that 
w <K t 1,dw , this probability is proportional to the accessible part of 
the illuminated area, which scales as w df . In particular, we expect 
that space- and time-dependence factorise, 

n(t;w)= ( d d rPUr,t) ~w di n (t), (8) 

J r<w 

where ITo(?) denotes the return probability to an infinitesimal vicin- 
ity of the origin. By the scaling law Eq. 0, we require that 
P m (x <sc 1) ~ 3 d ', which is confirmed by our simulations for 
the two-dimensional Lorentz model. As a by-product, one obtains 
n (0 ~ r d ' /dw = r d ' 12 , where d s = 2rf f /d w is the spectral dimen- 
sion. Combining both results, Yl{t\ w) ~ (wr 1/rfw ) rfr for sufficiently 
long times. 

Models of anomalous transport 

The Lorentz model. Anomalous transport emerges non-trivially 
in the Lorentz model.— ^ Here, a two-dimensional variant is used 
which consists of Brownian tracer particles exploring a disordered 
environment of randomly placed, overlapping circular obstacles of 
radius or, which we choose as a = 3 nm. The void space between 
the discs undergoes a continuum percolation transition at the criti- 
cal obstacle density n c cr 2 « 0.35907 (Ref. liol) . The infinite cluster 
displays self-similar behaviour characterised by the fractal dimen- 
sion d f = 91/48, known from lattice percolation. 25 The tracer dy- 
namics on this incipient infinite cluster is found to exhibit subdif- 
fusion, dri(t) ~ t 2/d ", with walk dimension d w ra 2.878 (Ref. [HI), 
see inset of Fig. [T] 

We have generated 1,600 trajectories of Brownian tracers with 
short-time diffusion coefficient Do = 2.5 um 2 /s, moving on the 



infinite cluster at criticality. (In practice, we computed trajecto- 
ries for particles on all clusters and evaluated the time-averaged 
mean-square displacement for each particle. Then, we selected 
those particles which did not show localisation based on a crite- 
rion for the local exponent of the mean-square discplacement at 
very long times; only these particles contributed to the final av- 
erage over independent trajectories.) Taking the divergent length 
scale into account, we have considered large systems of box length 
L = 10 4 cr = 30|im and have run the trajectories up to times of 
t ^ 10*%, where t = cr 2 /D = 3.6 (is is the natural time scale 
above which the diffusive motion is hindered by obstacles. The re- 
sulting correlation functions are invariant under time shift and do 
not display aging, in agreement with recent FCS experiments on 
crowded fluids. 1 For the in-silico experiment, we have evaluated 
the average in Eq. l(6j for beam waists between 24 nm and 384 nm. 

Fractional Brownian motion. Fractional Brownian motion 
(FBM) is a mathematical generalisation of the usual Brownian mo- 
tion yielding a subdiffusive mean-square displacement, dr 2 ^) = 
2dD a t a , with the generalised diffusion constant D a \ the distribution 
of the displacements AR(f) remains Gaussian. The description of 
a microscopic process generating such a dynamics is challenging, 
one formulation involving fractional derivatives was given in terms 
of a generalised Langevin equation. 32 Nevertheless, its "propaga- 
tor" (van Hove function) can be calculated exactly to 

Pimuir, t) = r" ? Gauss (re" 12 / ^ (9) 

where .PoaussW = (2^) _rf/2 x rf exp(-x 2 /2) and d denotes the dimen- 
sion of space. In particular, it satisfies the scaling form in Eq. 10 
exactly. Brownian motion with normal diffusion is obtained in the 
limit a — > 1, where D a becomes the diffusion constant. The FCS 
function corresponding to FBM is given exactly by Eq. ([5j. For 
comparison with the Lorentz model, we have fixed a and D a such 
that the mean-square displacements of both models coincide. 

Results and Discussion 

In the following, we will describe how FCS experiments with vari- 
able beam waist can provide insight into the microscopic dynam- 
ics and reveal spatially non-Gaussian, subdiffusive behaviour. We 
apply the generalised FCS theory from above to the exactly solv- 
able FBM model and to the two-dimensional Lorentz model with 
Brownian tracers. We have generated FCS correlation functions as 
described in the previous section. The obtained curves are shown 
in Fig. [T] and exhibit a significantly stretched decay compared to 
normal diffusion. For the corresponding FBM model with identi- 
cal mean-square displacement, the same trend, but a different shape 
of G(t) is found. For both models, an increase of the beam waist 
w shifts the relaxation to later times, while the shape appears to be 
preserved. 

Generalising the diffusion time in Eq. ([T}, we introduce the half- 
value time T i/2 (w) as a function of the beam waist via the implicit 
definition NG(t 1/2 ) = 1/2. The FCS data suggest a phenomenolog- 
ical scaling property, NG(t;w) = G(//ri/2(w)), i.e., all curves can 
be collapsed by appropriate rescaling of time. In the following, we 
will rigorously derive the scaling form of the FCS function G(t; w) 
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Fig. 2a: Simulated FCS data on double logarithmic scales; the data cor- 
respond to Fig. [T] The subdiffusive motion is mainly hidden in the tails of 
Goo(r; w) at long times, t/tQ 3> (w7<x) rfw , which decay as a power law with 
exponent a = df/d w , indicated by the broken line. 



Fig. 3a: FCS functions Goo(t) for various beam waists rescaled according 
to the scaling hypothesis, Eqs. (TTJ and )13t . Solid lines correspond to the 
simulation data of Fig. ^ dotted lines to the solution for FBM. Inset: The 
assumption of Gaussian transport, Eq. |5J, yields data collapse as well, but 
no saturation for large rescaled times. 
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Fig. 2b: Simulated FCS data as a function of the beam waist w for dif- 
ferent, but fixed time lags f; the data correspond to Fig. [T] The regime 
1 <*: w/<r <*: (f/fo)''"' w is characterised by a power law increase with expo- 
nent df (solid lines), revealing the fractal spatial structure. 

for the models under consideration. In particular, a thorough scal- 
ing analysis can discriminate whether or not a proposed theoretical 
model describes the spatio-temporal tracer dynamics contained in 
the FCS data. 

Scaling of the FCS function. 

The scaling properties of the FCS function are inherited from the 
van Hove function by means of the master formula, Eq. l|6}, 

NGUf, w) = J d d r exp (-r 2 /w 2 ) P m (r, t) . (10) 

In case of the Lorentz model, one finds from Eq. (|7) that 

NG m {t\w) = G m (wr Ud -) (11) 

for w » cr and t » to, and similarly for the FBM model. For both 
models, these scaling forms imply a power-law divergence of the 



Fig. 3b: The assumption of spatially Gaussian transport, Eq. (3), does not 
lead to scaling of the FCS data for obstructed motion. The anomalous ex- 
ponent a = 0.659 is obtained from the long-time decay of Goo(0 following 
Eq. {2). Dotted line: solution for the FBM model rescaled by the same 
procedure. 

FCS half-value time in particular, 

Ti/2(w) ~ w rfw for w»cr, (12) 

which is corroborated by the rescaling of G(t; w) below. For normal 
diffusion on a mesh grid model, a corresponding relation has been 
derived. — In the regime w <K t 1,dw , the FCS experiment essentially 
probes the return probability Tl(t;w). Thus G m (x « 1) ~ x*, and 
a non-trivial power-law decay of the FCS function is predicted at 
long times, 

G m (t;w) ~ r dfld ™ . (13) 

A double-logarithmic representation of our simulated FCS data in- 
deed renders the final decay of G„(?; w) straight lines, see Fig. [2a] 
Different beam waists yield parallel lines, and the slopes are com- 
patible with the expected value of df/d w = 0.66. 

The form of G„(-) for small arguments further yields the fractal 
space dimension, NG„(t\ w) ~ w' 1 ' for sufficiently large, fixed time 
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lag t. Thus, the structural properties can be obtained from the FCS 
data by fixation of the correlation time and sufficient variation of 
the beam waist, and the fractal dimension is directly accessible 
if the data are plotted on double-logarithmic scales, see Fig. [2b] 
The asymptotic regime is limited from below by the microscopic 
details of the system and from above by the crossover to the trivial 
behaviour, NG^it'^v — » oo) = 1. Within the remaining window, 
1 <k w/cr <sc (t/to) lldw , our FCS data reveal non-trivial power-law 
behaviour over 1.5 decades in space for the longest time lags, and 
the slope of the data corresponds to the fractal dimension d f = 1.9 
of the underlying space. 

A simultaneous test of both the temporal decay of G(t; w) and its 
dependence on the beam waist is provided by appropriate rescal- 
ing of the data. Figure [3a] shows the FCS function divided by 
the predicted long-time decay as a function of the rescaled time 
f = (t I t )(w I a-y d " . The excellent data collapse strongly corrobo- 
rates the scaling law of Eq. i ll It and of the half-value time, Eq. 4121 . 
Simultaneously, the rectified data converge for t — » oo, validating 
the algebraic decay, Eq. 4 13b . 

This analysis is in stark contrast with the form of G(t) obtained 
for subdiffusive motion and the assumption of spatially Gaussian 
transport, e.g., for the FBM model. It is instructive to discuss the 
implications for general dimension d. Then, Eq. {5} generalises to 

NGZ auss (t; w) = [l + Sri(t)lw 2 \'" 1 ~ w d r d/d ™ (14) 

as t — > oo, using Sr^if) ~ f 2/rfw . The inset of Fig. [3a] reveals that 
rectification with f 2/rfw does not lead to saturation at long times, re- 
flecting the fact that the decay of the FCS autocorrelation function 
for the obstructed motion decays with a different exponent than the 
inverse mean-square displacement. Nevertheless, the data still col- 
lapse on a single master curve, yet this shows merely that both the 
Gaussian ansatz and the critical scaling use the same reduced time 
i. Conversely for the FBM result, Gaussian scaling yields conver- 
gence at long times, while the assumption of obstructed motion 
does not. 

Let us perform the standard analysis for anomalous transport on 
our simulated FCS data, as suggested by Eq. l[2j. Fitting the long- 
time decay fixes the exponent to a = df/d w , which indeed yields a 
saturation in the rectification plot, see Fig. [3b] Interpreting this a 
as characteristic exponent of subdiffusion for the mean-square dis- 
placement, Sr^if) ~ f, and assuming Gaussian transport, Eq. 10, 
suggests the use of t ~ tw~ 21 " as scaling variable. However, then 
the data for our model system fan out for large times. Furthermore, 
it implies T\ii(yv) ~ w 2la = w 2dw,d, \ contradicting Eq. H2\ . We 
conclude that the standard approach using only a single scaling ex- 
ponent is not consistent for anomalous transport due to obstacles. 

We close with the question when the widely used Eq. <(2j is a 
valid description of the FCS correlation function. The equation is 
a specialisation of Eq. (0, which holds if and only if the distri- 
bution of displacements is Gaussian, i.e., solely determined by its 
second cumulant, 6r 2 (t)\ this is a consequence of the master for- 
mula, Eq. l[6). For spatially Gaussian transport, validity of Eq. {2} 
is then equivalent to a power-law increase of the mean-square dis- 
placement, 6r 2 (i) ~ f . If the functional form of Sr 2 (i) is different, 
e.g, if the dynamics exhibits a crossover from anomalous to nor- 
mal diffusion at some crossover time scale t x , Eq. $2% applies only 



to shorter time lags, t <K t x . Since such a crossover is generi- 
cally expected away from a critical point (see, e.g., the discussion 
in Ref. 18 and the simulation results in Refs. Ea - Jiij and liH) . the 
analysis of FCS data based on Eq. l[5) appears more robust. Finally 
as a test of the Gaussian assumption, it would be essential to quan- 
tify the corrections to Eq. lf5J by FCS experiments with variable 
beam waist. 

Conclusions 

We have shown that by systematic variation of the beam waist 
in FCS experiments, spatio-temporal information on the single- 
particle dynamics of complex systems can be collected. We have 
generalised the FCS theory beyond the assumption of spatially 
Gaussian transport and have derived a fundamental expression for 
the FCS correlation function, Eq. l[6j, which is a general starting 
point for the interpretation of experiments and which significantly 
facilitates theoretical and numerical work on FCS. In particular, it 
is straightforward to transfer our findings to the study of complex 
transport in other fields where FCS is widely employed, e.g., in 
physical chemistry and in polymer physics. 

The obtained master formula for FCS reveals an analogy be- 
tween FCS and time-resolved scattering techniques. It can be ex- 
tended to the case where the concentration of the labeled particles 
is not dilute any more. Then a distinct part arises in addition to 
the self-part similar to the corresponding decomposition of the co- 
herent intermediate scattering function. Likewise, one can easily 
account for the asphericity of the illuminated volume; yet this does 
not affect the scaling arguments presented in this work. 

For subdiffusive motion due to obstacles, both the fractal na- 
ture of the underlying space and the anomalous transport can be 
revealed by FCS. We have developed a scaling theory for G„(t\ w), 
which excellently describes our simulated data for the full range of 
investigated beam waists. These findings have been contrasted to 
fractional Brownian motion (FBM), an exactly solvable model for 
subdiffusion with different predictions for the scaling behaviour. 
The derived scaling properties should be experimentally accessi- 
ble with modern nanoscopic optical methods. 1314 In particular, the 
spatial information provided by FCS can be used to experimentally 
distinguish different routes to anomalous transport. 

We have demonstrated that a fit of the time-dependence of the 
FCS function for a single beam waist does not necessarily deter- 
mine the walk dimension, which characterises the subdiffusive in- 
crease of the mean-square displacement, 6r 2 (t) ~ t 2,d " . A more ro- 
bust procedure would be based on measurements of the half-value 
time Ti/2(w) of the normalised FCS correlations for a wide range 
of beam waists. Then, the exponent of subdiffusion may be ob- 
tained from the w-dependence Ti/ 2 (w) ~ w rfw , which is expected to 
hold for various models of anomalous transport, see Eq. dl2t . The 
knowledge of d„ would be the starting point for a scaling test of 
the full curves G(t ; w) similarly to Fig. [3a] which then would be 
completed by a characterisation of the decay at long times. The 
analysis by Wawrezinieck et al— is somewhat different from the 
one suggested here. Their "apparent diffusion time" T d f p corre- 
sponds to Ti/2 in our notation. Based on experimental data on the 
plasma membrane of COS-7 cells and on simulations, they find 
the phenomenological relation T1/2 ~ w 2 + const and discuss im- 
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Fig. 4: Experimental FCS half-value times ti/2(w) from FCS measure- 
ments on the transmembrane protein GfR-GFP of COS-7 cells in a double- 
logarithmic representation; the data were taken from Fig. 4D of Ref. [Tot 
The straight line indicates a power-law increase of the half-value times as 
function of beam waist, T1/2 ~ w 2 6 . 

plications of the offset as the data are extrapolated towards small 
beam waists, w — » 0. In Fig. [4] we have replotted their Ti/2-data 
for the transmembrane protein GfR-GFP as function of the beam 
waist on double-logarithmic scales. The data points nicely follow a 
straight line in agreement with our power-law prediction, Eq. d!2t . 
Considering the error bars and the limited w-range, we obtain an 
estimate of the walk dimension d w between 2.4 and 3.4; the most 
likely value is d w w 2.6, implying an exponent of subdiffusion of 
2/d w ss 0.77. We find this observation rather encouraging with re- 
spect to the applicability of our approach and to the usefulness of 
FCS with variable beam waist for addressing the leading questions 
on anomalous transport. 
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